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Abstract 



We study the "approximate squaring" map f(x) := x\x] and its behavior when iterated. We con- 
jecture that if / is repeatedly applied to a rational number r = l/d > 1 then eventually an integer 
will be reached. We prove this when d = 2, and provide evidence that it is true in general by giving 
an upper bound on the density of the "exceptional set" of numbers which fail to reach an integer. 
' We give similar results for a p-adic analogue of /, when the exceptional set is nonempty, and for 

r S^ iterating the "approximate multiplication" map f r (x) := r\x~\, where r is a fixed rational number. 

i > ■ 

C$ . We briefly discuss what happens when "ceiling" is replaced by "floor" in the definitions. 
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, In this paper we study the "approximate squaring" map / : Q — > Q given by 

CO 

P; f{x):=x\x], (1) 



1. Introduction 



and consider its behavior when iterated. Although there is an extensive literature on iterated maps 
(see for example [Collet and Eckmann 1980], [Beardon 1991], [Lagarias 1992]), including the study 
of various first-order recurrences involving the ceiling function ([Eisele and Hadeler 1990], [Graham 
J> ■ and Yan 1999]), the approximate squaring map seems not to have been treated before, and has 

^ , some interesting features. 

The function / behaves qualitatively like iterating the rational function R(x) = x 2 . Indeed, all 
points \x\ < 1 have a bounded orbit under f(x), while all points \x\ > 1 have unbounded orbits 
and diverge to oo, just as they do when R{x) is iterated. However, f(x) has the additional feature 
that it is discontinuous at integer points. It follows that the n-th iterate is discontinuous at a 
certain set of rational points, namely, those points x where f( n \x) is an integer. 

It is therefore natural to ask: if we start with a rational number r and iterate /, will we always 
eventually reach an integer? This question is the subject of our paper. 

Numerical experiments suggest that the answer to our question is "Yes" , although it may take 
many steps, and consequently involve some very large numbers. 

For example, starting at r = |, f(r) = | • 2 = 3, reaching an integer in one step; and starting at 
we get f(r) = ^, f^ 2 \r) = 41 > /^( r ) = taking three steps. On the other hand, starting 



r- 8 



7 

at r = |, we find 



,. . 12 m 36 , m . . 288 , f4V . 16704 , 5) 55808064 

fir) = — , f {2 \r) = — , fW{r) = — , f™ r = f { ) (r) = , 

5 5 5 5 5 
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, rfrt/ . 622908012647232 , 7 , , , 77602878444025201997703040704 

/ V '(»") = 5 . / l ( r ) = g , 

and we do not reach an integer until f^ 18 \r), which is a number with 57735 digits. 

We note that for any rational starting point r, since \x\ is an integer, the denominators dj 
of the iterates f^'(r) must form a nonincreasing sequence with dj+i dividing dj. For < r < 1, 
[r] = 1 and /(r) = r, so there the denominator is fixed. For —1 < r < 0, f(r) = 0, and for r < — 1, 
/(r) > 1. So it is sufficient to restrict our attention to the case of rationals r > 1. 

We make the following conjecture: 

Conjecture 1. For each rational r € Q mf/i r > 1, i/iere is an integer m > suc/i £/iat fv\r) is 
an integer for all j > m. 

We establish the conjecture in the special case when the denominator is 2, where a complete 
analysis is possible. This is done in Section 2. 

In Section 3 we consider the case of rational starting values r with a fixed denominator d > 3. 
We show that the set of starting values that reach an integer after exactly j steps has a limiting 
density, and that the set of starting values that never reach an integer has density zero. More 
precisely, at most a sparse subset of the rationals {r = k : d < I < x} fail to become integers under 
iteration, in the sense that the cardinality of this subset is bounded above by C(d,e)x 1 ~ ad+e for 
a certain positive constant ad and any positive e, where C(d,e) is a positive constant depending 
only on d and e. Showing that this "exceptional set" of starting values which fail to reach integers 
is in fact empty (or even finite) appears to be a difficult problem, for reasons indicated below. We 
also show that the set of starting values that reach an integer after exactly j steps has a limiting 
density. 

In Section 4 we consider a p-adic analogue of the approximate squaring map. In this case we 
show that there is a nonempty exceptional set of elements in which under iteration never 

"escape" to the smaller invariant set —^jWjp. This set has Hausdorff dimension exactly 1 — a p k, 
where a p k is the same constant that appeared in Section 3. The existence of this exceptional set is 
one reason why it may be a difficult problem to obtain better upper bounds on the cardinality of 
the exceptional set in Section 3. 

In Section 5 we study similar questions concerning the "approximate multiplication" map 

f r (x) :=r\x] , (2) 

where r is a fixed rational number. For r = t, this map is a special case of the map x i— ► a + [§] , 
where a,b £ Z, b > 2, studied by P. Eisele and R. P. Hadeler [Eisele and Hadeler 1990]. Recently, 
J. S. Tanton [Tanton 2002], together with Charles Adler, formulated a game-theoretic problem 
"Survivor" , and noted that its analysis leads to the study of the sequence of rational numbers 

a = r, a n = r \a n -{\ for n > 1 , 

for r > 1, which is the trajectory of r under the map f r (x). He raised the question, "Must some a n 
be an integer?", and conjectured that the answer is "Yes". This question differs from the case of 
the approximate squaring map in that the denominators of successive iterates, though bounded by 
the denominator of r, may increase or decrease. We note that the long-term dynamics of iterating 
this map differs according to whether \r\ < 1, |r| = 1 or |r| > 1, with the case r > 1 being most 
analogous to the approximate squaring map. 
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The approximate multiplication maps have some resemblance to the map occurring in the 3x + 1 
problem. Setting r = 4, we observe that f r (x) maps the domain ^Z into itself, and on this domain 
is conjugate to the map g r : Z — ► Z given by 



g r (n) 



if n = (mod d) , 
|(n + Z(d - 6)) if n = b (mod d),l<b<d-l, 



(see (44) of Section 5). In terms of the conjugated map the question we consider becomes whether 

3 
2 



for most starting values some iterate of g r is an integer divisible by d. For r = § the conjugated 



map is 



53/2 ( 



n 



|n if n = (mod 2) , 

3^ i 3 



if n + § if re = 1 (mod 2) 
This is similar in form to the 3x + 1 function 



T(n) 



in if 7i = (mod 2) , 

|n + i if n = 1 (mod 2) . 



as given in [Lagarias 1985] and [Wirsching 1998], although the long-term dynamics of g$/2 an d T 
are different. 

We formulate a conjecture for approximate multiplication maps analogous to the one above for 
the approximate squaring map. We define the exceptional set E(r) for the map f r to be 

E{r) := {n : n G Z and no iterate ff'(ri) £ Z for j > 1} . (3) 

Then we have: 

Conjecture 2. For each nonintegral rational r £ Q wii/i [r| > 1, i/te exceptional set E{r) for the 
approximate multiplication map f r is finite. 

The "expanding map" hypothesis \r\ > 1 is necessary in the statement of this conjecture, for 
the conclusion fails for r = ^ with d > 3, as remarked below. In parallel to the results for the 
approximate squaring map we prove Conjecture 2 for rational r having denominator 2; it remains 
open for all rationals with \r\ > 1 having denominator d > 3. This conjecture may also be difficult, 
as indicated in Section 5 by an analogy with the problem of showing there exist no Mahler Z- 
numbers, a notorious problem connected with powers of § ([Mahler 1968], [Flatto, 1991]; also 
[Choquet 1980], [Lagarias 1985]). 

Our analysis in Section 5 applies more generally to the family V T of maps h r : Z — ► Z having 
the form 

h r (n) = —(In + If,) when n = b (mod d) , (4) 
d 

where the integers If, satisfy = —lb (mod d). We obtain for all functions h r € V r an explicit 
upper bound on the cardinality of the exceptional set E(h r ,x) consisting of all integers \n\ < x 
that do not have some iterate that is a multiple of d. We show that for all rationals r, 

#E(h r ,x) < 4dx Pd , 

with (3d = lo fog^ • We complement this upper bound with a result showing that it is of the correct 
order of magnitude (to within a multiplicative constant) for certain values of r lying in < r < 1. 
This is the case for the function g r for r = -j with d > 3. It follows that Conjecture 2 cannot hold 
for these values of r. 
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Section 6 gives some numerical results related to these questions. 

As a final point we observe that one can also consider similar functions obtained by replacing 
the "ceiling" function with the "floor" function. Exactly the same set of questions can be asked for 
such functions. Our methods carry over to such functions, and there appear to be analogues of all 
our results and conjectures. For example: 

Conjecture 3. Let F(x) := For each rational r G Q with r > 2, there is an integer m > 

such that F^\r) is an integer for all j > m. 

We will not say any more about the "floor" function versions, apart from one result at the end 
of Section 5. 



Notation: ] denotes the ceiling function, [_ J the floor function, and {{j}- the fractional part. For 
a prime p, \ \ p denotes the p-adic valuation. (If r G Q , r = p a | with a, b, c G Z, c ^ 0, gcd(p, b) = 
gcd(p, c) = 1, then \r\ p = a.) Q p and Z p denote the p-adic rationals and integers, respectively. For 
integers r,s,i, r\s means r divides s, and r l \\s means divides s but r t+l does not. 



2. Denominator 2 



In this section we investigate the case when the starting value r has denominator 2. Here we are 
able to give a complete analysis. The following table shows what happens for the first few values of 
r. It gives the initial term, the number of steps to reach an integer, and the integer that is reached. 

f f .357 9111315 17 19 

Sldi l . 2 2 2 2 2 2 2 2 2 

steps : 1 2 1 3 121 4 1 • • • 

reaches : 3 60 14 268065 33 2093 60 1204154941925628 95 

The number of steps appears to match sequence A001511 in [Sloane 1995-2003] (and the num- 
bers reached now form sequence A081853 in that database). Indeed, we have: 

Theorem 1. Let r = with I > 1. Then f^ m \r) reaches an integer for the first time when 

m = \l\ 2 + 1. 

Proof. Note that if x G Q has denominator 2 and is not an integer then \x] = x + ^ • 

We use induction on \l\ 2 = v. If v = then / is odd, \r~\ = r + ^ = I + 1 is even, and r\r~\ has 

become an integer in one step, as claimed. 
Suppose v > 1, and 

I = 2 V + l v+1 2 v+1 + l v+2 2 v+2 + ■■■ 
is the binary expansion of /, where each Zj = or 1. Then 

r\r] = (z + iVz + 1) 

= I + 1 + / + /2 
2 2 

= \ + 2"" 1 + (Z„ +1 + 1)2* + (l v+1 + Z„ +2 )2*' +1 + • • • + 2 2v + ■■■ 
21' + 1 
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where 

l> = 2 V - 1 + (l v+1 + 1)2" + (l v+1 + l v+2 )2 v+l + • • • + 2 2v + • • • . 

and \l'\2 = v — 1. By the induction hypothesis, this will reach an integer in v — 1 steps, so we are 
done. ■ 

Remark. The numbers 21 + 1 for which \l\ 2 = v are precisely the numbers that are congruent to 
2 v+i + i ( mod 2 V + 2 ). For example, if v = 0, 21 + 1 G {3, 7, 11, 15, . . .}, of the form 3 (mod 4); if 
v = 1, 21 + 1 G {5, 13, 21, 29, . . .}, of the form 5 (mod 8); and so on. 

Corollary 1. Let r = ^f-,l > 1, \l\ 2 = v. Then the first integer value taken by f^ m '(r) is 



^ t,+1 )(2/ + l) , 



where 9(y) = y(y + l)/2. 



Proof. This is now a straightforward calculation, again using the fact that if x G Q \ Z has 
denominator 2 then \x~\ = x + \ . ■ 

For example, if v = 0, and r = (4A; + 3)/2 = y/2 (say), then in one step we reach the integer 
\0{y) = y(y + l)/4. If v = 1, and r = (8k + 5)/2 = y/2, then in two steps we reach the integer 

V%))-^ +1)(y2+y+2) - 



2 v 16 
if v = 2, and r = (16k + 9)/2 = y/2, then in three steps we reach the integer 

lo(e(9( y ))) = y(y + 1 Hv 2 + y + 2 ) (y 2 - v + 2 ) (y 2 + 3 v + 4 ) 

and so on. 

3. Denominator d 



We now analyze the case of rationals with a general denominator d, obtaining less complete 
results. The next theorem shows that most rationals will eventually reach an integer. In particular, 
it gives an upper bound on the number of such rationals below x that never reach an integer. Given 
an integer d > 2, and a bound x > 1 we study the "exceptional set" 

M d (x) := {I : l<l<x, f {m) ( 1 -) £ Z for each m > 1} , (5) 

and let M^(x) = |A4rf(a;)|. The finite set [1, d — 1] := {1,2, . . . ,d — 1} is contained in Md( x ), and 
Conjecture 1 asserts that M.d(x) = [1, d— 1]. 

Theorem 2. Let d > 2 be a fixed integer. 

(1) For each finite j > the set of positive integers I for which the approximate squaring map 
with initial value r = 4 first iterates to an integer after exactly j steps is a union of arithmetic 
progressions (mod cP +1 ) restricted to positive integers. 

(2) There is a positive exponent ad such that for each e > and all x > 1, 

M d (x) < C(d,e)x 1 ~ ad+£ (6) 
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for a positive constant C(d,e), with ad given by 



d'\d,d'>i \4>{d' 
where (j) is the Euler totient function. In fact 



a d = min log d , ( -— - ) , (7) 



a d = mm . (8) 

P j\\d jlogp 

Note: It follows immediately from (7) that < < 1 and ad = 1 only for d = 2. 

To prove the theorem we shall first prove a preliminary lemma. We need some notation con- 
cerning the pattern of denominators in a sequence of iterates. We write ^ = ^ with gcd(^o, do) = 1, 
and set 



d J dj 



where gcd(/j, dj) = 1 and dj\d with dj > 1. The pair (I, d) determines the sequence do, d±,d2, ■ ■ ■ ■ 

I 



We consider the first m iterates for a given input r = 4, and let 



y(l,d,m) := (do,di,...,d m ) (9) 

denote the vector of denominators of the first m iterates. We call this a chain of length m + 1. 
This information in 3^ can be described in another way. We set <i_i = d and define rj by 

r,- - 



for j = 0, 1, 2, . . . , m . We call j a break-point if > 1. Let the break-points up to the m-th iterate 
be < ji < j2 < • • • < jt < m. The data ^(Z, d, m) is completely determined by knowledge of the 
pairs (ii,^), (jt,rj t ) and m. Note that since each break-point removes a divisor of d from the 
denominator, in any chain of length m there are at most s break-points, where s is the number of 
prime factors of d, counted with multiplicity. 

Now consider an arbitrary chain y, which consists of a sequence (do, d%, d m ) with (ij+i di- 
viding dj and do dividing d. A chain is complete if d m = 1, and is incomplete otherwise. We will 
only consider complete chains with d m _\ > 2, i.e. with jt = in. 

The following lemma characterizes the set of positive integers / that have a given associated 
chain y of length m + 1. 

Lemma 1. Given a fixed chain y = (do, d\, d m ) of length m + 1, the set of positive integers I 
such that r = 4 has y(l, d, m) = y consists of a collection of arithmetic progressions 
(mod d-idodi ■ ■ ■ d m _i), restricted to the positive integers. There here are exactly 

^(d o )0(di)0(d 2 )"-0(dm) (10) 

such arithmetic progressions. 



(> 



Proof. We study the set of r = 4 having a fixed chain y for its initial m iterates. 

The set of allowable I with (I, d) = do consists of all integers r±lo with gcd(lo, do)=l. This 
consists of 4>{do) arithmetic progressions (mod d) = (mod ti-i). 

We now write ^- (for k = 0, 1, . . .) in a mixed-radix expansion where the radices depend on k: 

j- = + a (k) + £ Oj(*) II *W > C 11 ) 

ak ak j=l 1=0 

in which the "digits" a,j{k) satisfy 

< a_i(/c) < 4 , and < < dj + k for each j > . 

Here we set d m +j := d m for all j ' > 1. The sum on the right-hand side of (11) is actually a finite 
sum. By definition, 

h+i , ( h\ / a-i(k) 



dk+i 



^ (^) = ( ° 4^ + + ^ a i( k ) d kdk+i ■ ■ ■ 4+j-i 



1 + a (/c) + Y^a,j(k)d k d k+ i •• • 4+j-i J ■ (12) 

We use induction on k > to simultaneously establish four properties of this mixed-radix 
expansion: 

(i) We have 

r k+1 = gcd(a (fc) + 1,4) ■ (13) 



(ii) We have 



(iii) We have 



a_l(fe + l) = a„!(fc) ao(fc) + 1 (mod 4+i) • (14) 



gcd(a_i(A; + l),4+i) = 1 • (15) 

(iv) For < j < m, 

cij(k + l) = a-i(k)a,j + x(k) + G(a- 1 (k),a (k),a 1 (k), . . . ,a,j(k)) (mod d j+ k+i) , (16) 

for some function G(a_i(fc), ao(k), a\{k), . . . , aj(k)). This function includes all the necessary infor- 
mation about "carries" in the multiple-radix expansion. 

The base case k = is checked directly. Since 4 divides all terms in the sum on the right-hand 
side of (11), the right-hand side of (12) (when k = 0) has a single term a -i(°)(°o(o)+i) h av i n g a 

denominator, and this term equals a ~^ 1 ^ (mod 1). Since gcd(a_i(0), 4) — 1 by hypothesis, we 
must have 

gcd(a (0) + 1, 4) = -r = r i > 
«i 

which is (i) for k = 0. The term with a denominator in (12) is then 

ai{0) am±i 

4 
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hence 

a_i(l) = a_i(0) ( a ° ( ° r ) i +1 ) (mod dx) , (17) 

which is (ii). Now gcd(a_i(l), dx) divides gcd(a_x(0), d\) gcd( a °^ +1 , dx), both terms of which are 
1, so (iii) follows. Finally, to establish (iv) when k = 0, we drop the terms involving dxd2 • • • dj+x 
from (12) and observe that there is a term 

a_x(0)a J+1 (0)did 2 ---d j , (18) 

while all the other terms containing any o 3 -+i(0) for I > 1 are divisible by dxd 2 ■ ■ ■ dj+x- (Note that 
dxd 2 • • • dj + i divides d(jdxd 2 ■ ■ ■ dj.) This establishes that a congruence of the form (16) holds for 
the digit dj(l) in the expansion (12) for completing the proof of the base case. 

The induction step for general k follows using exactly the same reasoning. 

Next, we will deduce from (iv) that for each k > 0, 

k 

a (k + l) = a fc+ i(0)JJa_i(0 + G(o_i(0),ao(0),ai(0),...,o fc (0)) (moddjH-i) ■ ( 19 ) 
i=o 

for some function G depending on the indicated variables. The important point about (19) is 
that the dependence on afc+x(0) is linear, even though the dependence on the other initial terms 
a_x(0), ao(0), . . . , etfe(O) is nonlinear. 

To prove (19), we again use induction on k. For k = the assertion is that 

ao(l) = ai(0)a_i(0)+G(a_i(0),ao(0)) (mod d x ) , (20) 

which is (16) with k = j = 0. For k = 1 we wish to show 

ao(2) = O3(0)o_i(0)a_i(l) + G(o_i(0),ao(0),ai(0)) (mod d 2 ) . (21) 

Setting k = l,j = and fe = 0,j = 1 in (16) we obtain 

ao(2) = a_x(l)ai(l) + G(a_i(l),a (l)) (mod d 2 ) 

and 

ai(l) = a_ 1 (0)a 2 (0) + G(a_i(0),o (0),ai(0)) (mod d 2 ) , 

hence 

ao(2) == a 2 (0)a_x(0)a_i(l)+a_i(l)G(a_i(0),ao(0),ax(0)) + G(a_i(l),a (l)) (mod d 2 ) . 

However, from (17) and the fact that < a_x(l) < d\, a_x(l) is uniquely determined by a_i(0) 
and ao(0). Also the induction hypothesis allows us to use (20) to eliminate ao(l). Equation (21) 
follows. The case of general k follows in the same way; we leave the details to the reader. 
We have already seen that 

gcd(a_x(0),dx) = 1 . 

From (15) and (14) , 

gcd(a_i(2),d 2 ) = 1 , 
a_i(2) = a„x(l) Q ° (1) + 1 (mod d 2 ) . 

T2 
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Therefore gcd(a_i(l), cfe) = 1 and so 



gcd(a_i(0)a_i(l),d 2 ) = 1 • 

Continuing in this way we obtain 

k 

gcd([[a- 1 (!),d k+1 ) = 1 , (22) 

1=0 

for k = 0, 1, . . . , m. 

Then (19) shows that the congruence class of ao(k + 1) (mod dk+i) is uniquely determined by 
the congruence class of afc + i(0) (mod (4+i), once a_i(0), ao(0), . . . , flfc(O) are specified. 

In particular for each k > there are exactly (f>(dk+i) congruence classes of afc(0) (mod 4) 
that give 

,a (k) + 1 

gcd( = 1 , 

r k+i 

or in other words which give 

gcd(a (A:) + 1,4) = r k+i ■ 

At each iteration we impose one such condition. Combining this with the (j)(do) congruence classes 
(mod d-i) allowed initially, (10) follows. ■ 



Proof of Theorem 2. Part (1) follows from Lemma 1. Indeed the set of integers that iterate 
to an integer in exactly j steps are precisely those integers that belong to some complete chain 
y having j = jt. There are finitely many complete chains of length m + 1 with m = jt = j, 
and by Lemma 1 the set of positive I such that K belongs to a given such chain is a union of 
arithmetic progressions (mod d^\d^di . . . d m _i). Each such progression subdivides into a finite 
union of arithmetic progressions (mod d m+1 ), since d-id$d\ . . . d m _\ divides d m+l . 

To establish part (2), let the cutoff value x be given. We call a value I < x "bad" at size x if 
there is an iterate m > such that 

m— 1 m 

Yl dj < x < Yl dj , (23) 

3=0 3=0 

with d m > 1; here m depends on I, d and x. The chain y associated to the data (do, di, d m ) is 
necessarily an incomplete chain, and we say it is incomplete at size x if (23) holds. Conversely all 
elements I < x belonging to such an incomplete chain are "bad". We let N^x) denote the number 
of "bad" elements at size x, and we clearly have 

M d (x) < N d {x) . 

We will establish part (2) by showing the stronger result that 

N d (x) < C(d,e)x 1 ~ ad+£ . (24) 

To bound N d (x) we first bound the number of distinct incomplete chains at size x, and then 
bound the maximal number of elements below x falling in any such chain. It is easy to bound the 
number of incomplete chains y at size x as follows. Since each dj > 2 we have m < log 2 x. There 
are (™) choices for the break-points, and the bound t < s then yields that there are at most 

£( log t 2l )< s (iog 2 z)' 
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break-point patterns. Since a nontrivial divisor n of di-\ is removed at each break-point, we 
conclude that the number of distinct chains is at most d\s(log 2 X Y ■ F° r d (and hence s) fixed this 
gives for any fixed e > and x > 2 that 



#(y)<C 1 (d,e)x s . 



(25) 



for some positive constant C\(d,e). 

Now let N d (y,x) count the number of bad I at size x having a given chain y = (d$,di, ...,d m ) 
that is incomplete at size x. In view of (25) it suffices to prove an upper bound of the same form 
(24) for each such chain. So consider a fixed such chain y. Now (23) holds, and shows that each 
arithmetic progression (mod cfodi ■ ■ ■ d m -i) contains at most d m elements below x. From (10), 
Nd(y,x), the number of bad elements I at size x with chain y , satisfies 



Now x > dodi 



N d (y,x)<d m <i>(d 1 )c/)(d2) 

d m -i > did? ■ ■ ■ d m -i, hence 



4>(dr, 



N d (y,x) 



< 



dm J 



<t>{dj) 



dj 



and therefore 



N d (y,x) 



< d m e 



< 



dx~ 



log X 



Now x < dodi ■ ■ ■ d m , so that 

V m log- ^ 
Z^=o i0 § d 



<t>(dj) 

■j 



logx 



> 

> 

> 
> 



E 



m 
3=0 



log- 



mm 

0<j<m 



mm 

0<j<m 



log ^) 



log dj 



<f>(dj) 



where we used the definition (7). Substituting this in (26) yields 



(26) 



since , 3 < 1 gives log , < 0, and so 

cij <ij 

N d (y,x) < dx 1 -^ , 

as required. Combined with (25), this gives (24), hence (6). 
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Finally, we establish the equivalence of (7) and (8). Since 4> is multiplicative, for a general 
d! > 1 we have 



loe 



pi\\d> iu £> \4>{p>) 



> min (log^ ( J 1 . (27) 

Thus the minimum in (7) is attained when d' is a prime power. Now 



log 



Knp 1 )/ log jiogp 

is minimized by making j as large as possible, so we obtain the formula (8). ■ 

The upper bound in (24) in Theorem 2 has an exponent that is essentially best possible. Indeed, 
if l,d,x,m are such that 

d\ = d 2 = ■ ■ ■ = d m = d' , 

where d! is the value that minimizes (7), so that t = 1 and the chain is simply y = {(l,d')}, we 
have 

N d (y,x)>c' d x 1 - a ^ , 

for some positive constant C d . 

We next give a probabilistic interpretation of Theorem 2, following a suggestion of Daniel 
Berend. Given a denominator d and an integer cutoff value x, let X be a random integer chosen 
uniformly on the interval [l,a;], and let s = X/d. We iterate / starting at s, and let Stop(X) 
be the random variable giving the number of iterates needed to reach an integer. Stop(X) takes 
nonnegative integer values (including and oo). There is an associated probability distribution 

{Prob^Stop = j] : j = 0, 1, . . . , oo} 

which depends on the interval [l,x]. 



Theorem 3. For each fixed denominator d, there exists a limiting distribution as x — > oo of the 

I 



stopping times of the approximate squaring map f on rationals k. More precisely, for each finite j 



the limit 

Prob^StopiX) = j] := lim Prob x [Stop{X) = j] (29) 
exists and is a rational number with denominator dividing d J+ , and for j = oo we have 

Proboo[Stop(X) = oo] := lim Prob[Stop(x) = oo] = . (30) 

x— >oo 

In the case when d = p is a prime the limiting distribution is a geometric distribution: 

Prob^ [Stop(X) = j] = p(l - P y . (31) 



Proof. The assertion (29) follows directly from part (1) of Theorem 2, which shows that the set 
with stopping time j is a union of complete arithmetic progressions (mod d J+1 ). The assertion 
(30) follows from part (2) of Theorem 2. The computation of the density (31) follows from the fact 
that for d = p and each j > there is a unique complete chain y with stopping time j, having the 
single break-point (J, 1), giving one arithmetic progression (mod p ?+1 ). ■ 
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Remarks. (1) This stopping time formulation is similar to that arising in the 3x + 1 problem 
and studied in [Terras 1976, 1979]. 

(2) In principle the limiting distribution for any specific composite d is computable from the 
proof of Theorem 2. 

(3) For all d we have Prob 00 [Stop(X) = 0] = 1/d. 

4. p-Adic iteration 



We now consider an approximate squaring map defined on the p-adic numbers analogous to the 
approximate squaring map on Q, and study the question of whether some iterate will eventually 
become a p-adic integer. We show that now there is always a nonempty exceptional set of p-adic 
numbers which never become p-adic integers. 

Let p be a prime > 2 and let Q p denote the p-adic numbers, with typical element a = 
^2JL-k a jP J ^ where k € Z and the a,j satisfy < dj < p — 1. The p-adic integral part of a is 
given by the function F p : Q p — > Q p , where 

oo 

F p (a) := ^2 a iP j ' 

j'=max{0,fc} 

while the p-adic fractional part (or "principal part") of a is 

P p (a) :=^2a j p j , 

which is a finite sum (possibly empty); thus a = P p (a) + F p (a). We investigate the function 

f p (a) := a{F p {a) + 1) . 

which is ap-adic analogue of the approximate squaring map defined in (1). If we regard the rationals 
Q as embedded in Q p , then for nonintegral r in the subring Z[^] CQC Q p , we have f P (r) = r\r], 
so the iterates there agree with the approximate squaring map. 

For each k > the set Z p is invariant under the action of f p (x), and 

ZC-ZC \z C ••• . 
p p l 

We define the exceptional set 

fifc(p) := |a G ^Z p : f^ m \a) £ -L Z p for each m > l| , (32) 

for A; > 1. The set fifc(p) is an analogue of the exceptional sets Md(x) studied in the last section, 
corresponding to the denominator d = p k . (Note that Clk(p) H \%>o is contained in \M. p k(oc), as 
defined in (5)). We will show that these sets are nonempty, determine their Hausdorff dimension 
s, and get upper and lower bounds on their Hausdorff s-dimensional measure. 

The s-dimensional p-adic Hausdorff measure 7i p {Vl) of a closed set in Q p is defined by the 
general prescription in [Falconer 1990, Chapter 2] or [Federer 1969, Section §2.10]. Here < s < 1. 
The diameter of a measurable set S C Q p is 

\S\ := sup{|o — 0\ p : a, f3 G S} , 
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and its p-adic measure fi p (S) is Haar measure with the normalization fi p (Z p ) = 1. A p-adic ball 

B(a;p l ):={f3E® p : \(3 - a\ p < p 1 } 
is both closed and open, and has the property that its diameter equals its measure: 

\B{a ] p l )\=ix(B{a ] p l )) = - j 

For each 5 > we define 



p l 



H s {n, 5) := inf I ^ \Ij\ s : fic[J Ij , \Ij\ < 5 for all j \ , (33) 

and 1 

n s Jn) -.= HmH s (n,s) . (34) 

The p-adic Hausdorff dimension of fi is the unique value sq such that 7ip{0) = oo for s < sq and 
7i.p(ft) = for s > so- The value of 7^p°(0) may be zero, finite or infinite. 

A closed set f] C Z p is called weakly self-similar (mod p k ) with branching ratio b, where b is an 
integer > 2, if the following "equal branching" property holds for I = 1,2, . . .. Let Wi(£l) denote 
the set of initial sequences of digits of length Ik in Q, i.e. 

Ik-l 

Wi(9) = {P =^2 a i p ' : tliere exists some a G ^ whose "initial part" p lk P p (p~ lh a) = (3) . 

3=0 

(35) 

Then each sequence in Wi{Q) should extend to exactly b sequences in Wi + i(Q). That is, if the 
digits of an element (3 6 £1 are grouped in blocks of size k, once the first I blocks of digits are 
specified, there are exactly b allowable choices for the next block of digits. 

Theorem 4. For each k > 1 the set p k Q,^{p) is weakly self-similar (mod p k ) with branching ratio 
b = 4>(p h ) =p k - p fe_1 . 

Proof. The set Wi{£lk{p)) specifies the conditions under which the first I iterates /(a), f^ 2 \a), 
. . ., f( l \a) ^ Z p . Each such condition is a congruence (mod p lk ), which has exactly (fi(p k ) = 
pk _pk-i solutions for the next digit (compare (16) in the proof of Theorem 2, taking each dj = p k ). 
Finally the definition of Qk(p) implies it is a closed set. Thus p k Qk(p) is weakly self-similar. ■ 

We can determine the Hausdorff dimension of weakly self-similar sets in Z p , together with upper 
and lower bounds for the Hausdorff measure at this dimension. 

Theorem 5. Let Q, C Z p be a weakly self-similar set (mod p k ) with branching ratio b satisfying 
2<b<p k . Then Q. is a compact set and has Hausdorff dimension s(S7) given by 

logp K 

Its s(O)- dimensional Hausdorff measure satisfies 

h 

< n;^(n) < 1 . (36) 



1 The limit (which may be oo) exists since 7i s (fl, <5i) > H. a {Q,,&2) if <5i < S2 
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Proof. The compactness of is established similarly to Theorem 4. We define W/(Q) as in (35). 
This set has cardinality b l by hypothesis, and 

oo 

n = f]Wi(n) , 

1=1 

where Wi(£l) is the compact set 

W l (Q) = {p€Z p : = (3 (mod p lk ) for some (3 £ Wj(fi) | . 

To determine the Hausdorff dimension it suffices to establish the inequalities (36), since the 
fact that the Hausdorff measure is positive and finite determines the Hausdorff dimension. Set 
s = (log6)/(logp fc ), so that p ks = b. Also < s < 1. 

For the upper bound in (36) we consider the sets Wi(Q). Now b l balls of diameter p lk cover 
Wi(Q), hence cover J7. Thus for 5 = p~ lk this covering gives 

hi U 

which implies 

The lower bound argument is similar in spirit to that used for Cantor sets in [Falconer 1990, 
pp. 31-32]. By the compactness of we need only prove that the lower bound holds for finite 
coverings. The non-archimedean property of the valuation | \ p means that each Ij has diameter p m 
for some m, and hence we can enlarge Ij to a ball B{a;p m ) 2 Ij without changing its diameter. 
But B(a;p m ) gives an open cover, so it has a finite subcover: 

in 

3=1 



and 



We want to show 



m 

rrijS 



3=1 3=1 



k-l 

b \ k 



T,\ B (^p mj )\ s * (- k ) ■ ( 37 ) 

3=1 XP J 

We first replace these balls with balls of diameter p~ kl i for integers L. Write 

rrij = —klj + kj , < kj < k — 1 . 
Then we can cover B(aj;p mj ) with p k i balls B{aji\p~ kli \ We claim that 

m / h \ TT" 

^|fl(a i; j^)|' > K £|B(a,v;p-«0|". (38) 

3=1 KP J j' 

This will follow if we show that 
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holds for each j. This in turn is equivalent to showing 



Since p ks = b we have 



b \ k I b\ k 



\P J \P 



which proves (38). Thus (37) will follow from showing 

m 

Y,\B(a 3 ;p- kl i)\ s > 1 (39) 
i=i 

for any set of such balls that covers f2. We may suppose l± < li < ■ ■ ■ < l m . By the weak self- 
similarity of 17 there are b lm principal parts to cover with balls of diameter p~ klm in W m (Q). Weak 
self-similarity also says that each ball of radius p~ kl i covers either none or else exactly b lm ~ l i such 
principal parts in W m (fl). Since the balls cover $7, we must have 

m 



Dividing by b lm yields 



Using this bound we obtain 



b~ l i > 1 . (40) 

i=i 



m 



Y,\B{a f ,p-^)\ s = £p- 

3=1 3=1 



m 

u 



> 1 , 



which is (39). Thus (37) holds. ■ 

Corollary 2. The set Vtk{p) has Hausdorff dimension 



kg ( 1 + 



dim H (n k (p)) = s(p k ) := 1 ) (41) 

fclogp 



for k>l. Furthermore, its s(p k ) -dimensional Hausdorff measure satisfi 



es 



(i-^) 1 V-p* -1 ) < nf k \n k ( P )) < p k -p k ~ l . (42) 
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Proof. This follows by applying Theorems 4 and 5 to the set p k Q.k{p) and using the fact that 

H s p ^\n k (p)) = bH^ k \p k n k ( P )) , 

since (p k ) s = b. Using the branching ratio b = 4>(p k ) = p k — p k ~ l , we have 

logft logp k + log(l — 1/p) 
logpfc logp fc 

which gives (41). ■ 

Remark. For the branching ratio b = 4>(p k ) = P k ~ V > one can show that equality may occur 
on either side of the Hausdorff measure bounds in (36) (or (42)). For the lower bound, take the 
p k — p k ~ 1 allowed digit sets in each layer to be Yl^=ik 1 a j P' with the restriction that aik ^ 0. 
Then we can cover f2 with (p k — p k ~ l ) l+l (p — 1) balls of radius p~( lk+1 \ and get 

H s p (n k ( P )) < ( P k - p k ~ l ){ P k - p k ~ l ) l p- lks [(p - i)p~ s \ 

Since this coincides with the lower bound in (42), Hp(£lk(p)) must equal this bound. 

For the upper bound, choose the p k — p k ~ l allowed digits in each layer to be Y^j=tk * a i P* 
with the restriction that ak(i+i)-i 0- Then each residue class (mod p- k ( l + 1 )+' r "j covers exactly 
(p ~ l)p r_1 classes of (mod p~ lk ), and we can do no better than the upper bound. 

It would be interesting to know how many elements r = -j^ have no iterate /" (r) € Z p . We 

know that this set contains the p k — p k ~ 1 elements 1 < I < p k — 1 with (l,p) = 1. We conjecture 
that these are the only such elements. 

5. Approximate multiplication maps 




We can use similar methods to study iteration of the approximate multiplication map f r : Q — > Q 
given by 

fr(x)=r\x], (43) 

where r is a fixed rational number, say r = K with gcd(l,d) = 1. In this case, since x enters into 
the iteration only as the integer \x~\ , we may restrict attention to initial values i£Z. We consider 
the case when the denominator d > 1, and study the question of whether some iterate f^\x) will 
be an integer for some j > 1. Note that all iterates lie in ^Z. 

Unlike the case of approximate squaring, the iterates do not remain integral once they become 
integral. However, the truth of Conjecture 2 would imply that infinitely many members of a 
sequence of iterations {f^\n) : j > 1} will be integers, provided \r\ > 1. 

It is convenient to rescale the map to eliminate the denominators d, by conjugating f r by the 
dilation <&d(x) = dx. The result is the map g r ■ Z — > Z given by g r (x) := &d° fr ° ^ 1 (x)- Thus 

gr{x) = df r £) = l\*-\ , (44) 
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with 



for j = 1,2, We have 
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U) 



(x) = d/W(^) 



<7 r (n) = ^(^ n + ^) when n = b (mod d) , 

where Iq = and 

l h = l(d-b) for 1 < 6 < d - 1 . 
For example, when r = |, we have 



53/2O) 



|n if n = (mod 2) , 

\n + § if 77, = 1 (mod 2) . 



Our question then becomes: when does the sequence of iterations {gP\n) : j > 1} contain a term 
which is divisible by d? 

The map g r belongs a general class of functions which we will denote by V T , consisting of those 
"periodically linear" functions h r : Z — > Z of the form 

h r (n) = —(In + when n = b (mod <i) , (45) 
d 

where r = l/d is rational and the integers {lb : < b < d — 1} satisfy the conditions 

Z 6 = -Z6 (mod d) (46) 

needed to give an integer-valued map. Thus, although the notation does not reflect this, h r is 
defined by specifying r = l/d and constants lo,h, ■ . ■ , ld-i satisfying (46). We note that any map 
h r in V r is "self-similar" in the sense that its linear part is independent of the residue class. 
More general classes of periodically linear functions have been studied in connection with the 3x + 1 
problem - see §3.2 of [Lagarias 1985]; those classes differ from V r in allowing the linear part of the 
map to depend on the residue class (mod d). 

Our methods apply generally to the question of whether a particular function h r € V r has 
an iterate that is divisible by d (the denominator of r). The behavior of the map h r depends 
on whether |r| > 1, the "expanding map" case; r = ±1, the "indifferent map" case; or |r| < 1, 
the "contracting map" case. Our motivation comes from the expanding map case, but the results 
proved below apply to all cases. 

We formulate a general conjecture concerning functions in this class that are expanding maps, 
of which Conjecture 2 is a special case. 

Conjecture 4. Let r = ^ with gcd(/, d) = 1 and \r\ > 1. Let h r : Z — > Z be a function in the class 
V r . Then for each integer n, with at most a finite number of exceptions, there is some iterate j > 1 
such that hr(n) = (mod d). 

The "expanding" condition on r is necessary, for the conjecture fails for certain r with < r < 1, 
as shown in Theorem 8 below. 

An interesting map in the class V r is 

g r (x) = \rx\ . (47) 
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This map has Iq = and 

l b = d- {lb mod d) for 1 < 6 < d - 1 . 

The function 93/2(2;) appears in Mahler's study of Z-numbers [Mahler 1968], as explained below. 
We have 



53/2 (n) 



\n if n = (mod 2) , 

|n + i if n = 1 (mod 2) . 



Mahler's study of Z-numbers led to questions similar to Conjecture 4. A Z-number is a positive 
real number £ with the property that 

< t(|m < ^all n>l, 

where |xjf = x — [xj denotes the fractional part of x. Mahler conjectured that Z-numbers do not 
exist, and showed that a necessary and sufficient condition for their non-existence is that for each 
n > 1 there exists some j > 1 (depending on n) such that 

~g^ 2 {n) = 3 (mod 4) . (48) 

Mahler's conjecture remains open. Mahler obtained a nontrivial upper bound on the number of 
Z-numbers smaller than x, and [Flatto 1991] improved the upper bound to 0(x 0,59 ) for x — > 00. 

In comparison to (48), Conjecture 4 for g 3 /2 (with r = |) asserts that for each n £ Z there 
exists some j > 1 (depending on n) with 

g { 3 %(n) = (mod 2) , (49) 

with at most a finite number of exceptions. This special case of Conjecture 4 is true, as a conse- 
quence of the next theorem. There is exactly one exceptional integer, n = — 1, whose iterates never 
satisfy (49). 

More generally, for the case of rational numbers r with denominator d = 2, Conjecture 4 is 
provable for all functions in the class V T in a fashion analogous to that used for the approximate 
squaring map in Section 2. 

Theorem 6. Let h r be a function in the class V r , for fixed r = where t is an integer. Then 
for each n 6Z, with at most two exceptions, there exists some iterate k > 1 with 

hi k) {n) =0 (mod 2) . (50) 

Proof. We have 



h r (n) 



rn + lo if n = (mod 2) , 

rn + ^ + Ii if n = 1 (mod 2) . 



for some integers Zo, 

We claim that the set of integers n satisfying 

h®(n) ^ (mod 2) for 1 < j < k (51) 

consist of the integers in exactly two arithmetic progressions b (mod 2 fc+1 ), one consisting of even 
integers and one of odd integers. We prove the claim by induction on k > 1. For the base case, let 
b = oq (mod 2) with a$ = or 1 fixed, and consider the arithmetic progression n = b + 2m, with 
m £ Z. Then 

h r {n) = h r (b) + (2t + l)m , 
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and the condition h r (n) = (mod 2) then restricts m to lie in a single congruence class m = 
a\(b) (mod 2). We conclude that exactly two congruence classes b = ao + 2a± (mod 4) satisfy (51) 
for k = 1, one consisting of even integers and one of odd integers, completing the base case. 

For the induction step, supposing (51) true for k, let b (mod 2 fc+1 ) run over the two allowed 
congruence classes for the given k. Consider the arithmetic progression n = b + 2 k+1 m, with m € Z. 
Then we have 

h^\n) = hl j \b) (mod 2) for 1 < j < k , 

and 

4 k+1 \n) = h k r +1 (b) + (2t + l) k+l m . 

The condition that 

h\ k+l) (n) = (mod 2) 
is equivalent to m = ak + i(b) (mod 2), which excludes the congruence class 

b' = b + a k+l (b)2 k+l (mod 2 k+2 ) . 

Thus two congruence classes (mod 2 k+2 ) remain which satisfy (51) for 1 < j ' < k + 1. Since each of 
the previous classes b (mod 2 k+1 ) contributed one of these classes, one contains even integers and 
the other contains odd integers. This completes the induction step. 

Denote these two classes by bo(k + 1) (mod 2 k+2 ) and b\(k + 1) (mod 2 k+2 ), respectively. It 
follows that (50) holds except for integers in the sets 

oo 

\J{n = b {k) (mod 2 fc+1 )} . 

k=l 

and 

oo 

(J{n = b ± (k) (mod 2 k+1 )} . 

k=l 

Each of these sets contains at most one element, so there are at most two exceptional elements. ■ 
Concerning Theorem 6, there exist examples of functions h r G V T with denominator d = 2 such 
that there are zero, one or two elements in the exceptional set. For example, for r = | the function 
<7 3 / 2 has the single exceptional element n = —1. The map 

u i ^ / l n + 5 if rc = l(mod2), / K<Vl 
h ^ n) = \ln-l ifn = 0(mod2) (52) 

has two exceptional points, and —1, all the iterates of which are odd. 

We also note that Theorem 6 applies to the two values r = ±^ where the maps in V T are 
contracting. 

We next prove a result in the direction of Conjecture 4 for denominators d > 3. We bound the 
number of exceptions below x for a general function in the class V r , using an argument similar to 
that for the approximate squaring map studied in Section 3. Given a rational number r = 4 with 
gcd(Z, d) = 1 and a function h r in the class V r , we define the exceptional set by 

E(h r ;x) :={neZ: |n| < x, h^ k \n) ^ (mod d) for all k > 1} , (53) 

and let 

N(h r ;x) :=#(E(h r ;x)) . 

(The exceptional set in (3) is E(r) = E(g r ;oo).) The following result holds for all rational r, 
including those with |r| < 1. 
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Theorem 7. Let r = 4 be a rational number with gcd(7, d) = 1, and suppose that d > 2. 

^ .For eac/i function h r G i/ie sei o/ integers I such that the function h r iterated starting 
from the initial value x = I reaches an integer for the first time at iterate j > 1 is a union of 
d(d — iy complete arithmetic progressions (mod d^ +l ). 

(2) There is a constant < (3d < 1 depending only on d such that for every function h r G V r 
we have 

N{h r ;x) < 4dx A . (54) 

The precise value of the constant is 

^^--'°^ + ^>- <«> 

Proof. The arguments to establish this result are simpler than those for the approximate squaring 
map because we can use radix expansions to the fixed base d. If n = b (mod d) then 

h r (n) = l -n + ^ , (56) 

where < lb < d — 1 with lb = —lb (mod d). 
We claim that the elements n G Z such that 

h®(n) ^ (mod d) for 1 < j < k (57) 

consist of a certain set of d(d — l) k residue classes (mod d k+1 ). We proceed by induction on k > 1. 
For the base case k = 1, given 6 (mod d) the elements of the arithmetic progression n = b + dm 
with m e Z have 

h r (n) = h r (b) + /m . 

Since gcd(/,m) = 1, as m G Z varies these numbers cycle through every residue class (mod d). 
In particular there is one class m = ai(b) (mod d), say, that gives g r (n) = (mod d). This 
arithmetic progression b + a\{b)d (mod d 2 ) is ruled out and all elements of the remaining d — 1 
arithmetic progressions (mod d 2 ) satisfy (57) with k = 1. This completes the b asc case. 

For the induction step, suppose (57) holds for k, and there are d(d — l) k allowed residue classes 
b (mod d k+1 ). For each of these residue classes consider the arithmetic progression n = b + d k+1 m 
with m € Z. Using (56) repeatedly, we have 

h k+1 (n) = 4 k+1 \b) + l k+1 m . 

Since gcd(/,m) = 1 this progression cycles through all residue classes (mod d), and the condition 
h k+1 (n) = (mod d) rules out one residue class b + af c +i(b)d k+1 (mod d k+2 ), say. Thus imposing 
the additional condition 

h k+1 {n) £ (mod d) 

leaves d(d — l) k+l allowed residue classes (mod d k+2 ) whose elements satisfy (57) for k + 1. This 
completes the induction step, proving the claim. 

The proof of the claim establishes part (1), since the elements that first have an integer iterate 
at the k-th step for k > 1 were shown to form d{d — l) k ~ l complete arithmetic progressions 
(mod d k+l ). 

To establish part (2), we will prove a stronger result than stated above. Set 

E*{h r ;x) := {n G Z : \n\ < x, h {k \n) ^ (mod d) for 1 < k < log d x} (58) 
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and let N*(h r ;x) = #(E*(h r ;x)). Certainly N*(h r ;x) > N(h r ;x), so it suffices to establish 

N*{h r -x) < Adx Pd , (59) 

where (3^ is given in (55). 

Now suppose d k < x < d k+1 . We observe that applying (57), with k replaced by k — 1, shows 
that the elements in E*(h r ;x) are necessarily contained in a set of d{d — l) fc_1 residue classes 
(mod d k ), since log d x > k — 1. It follows that 

N*(h r ;x) < d(d- I)*" 1 -2 

since <i > 2. Now x > d k , so k < and we have 

N*(h r ;x) < Ade kXo ^ d -^ 

log(rf-l) 

as asserted. ■ 

Theorem 7 has a probabilistic interpretation analogous to that of Theorem 3. We leave its 
formulation to the reader. 

We now show that the upper bound of part (2) of Theorem 7 is nearly best possible for certain 
values of r in the interval < r < 1, where the map is a contracting map. That is, for suitable 
maps in the class V r , the upper bound (54) of Theorem 7 is within a multiplicative constant of the 
best possible upper bound. 

Theorem 8. Let r = k with d > 3. Then, for the conjugated approximate multiplication map g r 
of (44) > the exceptional set E(g r ;x) has cardinality 

N(g r ] x) > -x^ d , for all x > d , (60) 
with Pd = 10 Tog 1 ^ • ^ n V ar ticular, the full exceptional set E(g r ; oo) is infinite. 

Proof. Note that for r = -j the functions g r and g r coincide. We claim that, for each k > 1, the 
subset Sfc of [l,d k ] given by 

{n : 1 < n < d k , n = oq + aid + . . . + afc_i<i fc_1 with all ^ —1 (mod d) for i > 1} 

is contained in E(g r ;d k ). (The cardinality of £/% is d(d — l) fc_1 .) We prove this by induction on 
k > 1. 

For the base case A; = 1 we have 

Si = {l,2,...,d-l}C^ r ;d) 
because each <7 r (n) = 1^0 (mod d) and 1 is a fixed point of g r . Next, 

5r(Sfc) C , 

because 

g r (n) = {a\ + 1) + a 2 d + a 3 d 2 H h ak-\d k ~ 2 G Sfe-i , 



< 2<i 2 ((i — l) fc_1 < 4d(d- 1) 
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and l<ai + l<d— lby hypothesis. Thus C E(g r ; d k ), which completes the induction step. 
For d k < x < d k+l , with k > 1, we have 

N(g r ; x) > N(g r ; d k ) > #S fc = d(d - l)*" 1 > -x^ , (61) 

d 

as asserted. ■ 

Since the exceptional set E(r) of the approximate multiplication map f r has E(r) = E(g r ; oo) 
and E(g r ,x) C E(g r ;oo), Theorem 8 shows that the conclusion of Conjecture 2 does not hold for 
these values of r. 

It seems plausible that for all values —1 < r < 1 (except r = 0) there is some function in the 
class V r for which the conclusion of Conjecture 4 does not hold. However we do not attempt to 
construct such functions here. 

We end this section with a result that shows that in at least one case the approximate multi- 
plication map based on the floor function, 

F r {x) := r[x\ , 

behaves in exactly the same way as f r (x). A version of this result was communicated to us by 
Benoit Cloitre. 

Theorem 9. Let r = with d > 1. Then, for any integer m > 1, the sequence iy^(m + d), 
j > 0, takes exactly as long to reach an integer as as the sequence 

Proof. Let y i = $\m) and Yj = F^\m + d). An easy induction argument shows that Yj = 
yj + d + 1 for j > 0. ■ 



6. Numerical results 



The simplest case where we do not know if the approximate squaring map / of (1) will always 
reach an integer is when the starting value r = l/d has denominator d = 3. We wish to determine 
6{r) (say), the smallest value of k > for which f^ k \r) is an integer. 

Testing any particular value of r is complicated by the fact-already illustrated in Section 1— 
that the iterates grow so rapidly. This difficulty can be overcome by writing the k-th. iterate 
h/d k ■= f [k) {r) in "base d": 

, oo 

f = £ a 3 {k)d k , (62) 

where the "digits" dj(k) satisfy < aj(k) < d (compare (11)), but storing only the terms in (62) 
with j < M. That is, we work mod d M+1 . As long as 6{r) < M — 1, we get the correct answer by 
finding the smallest k for which a~i(k) = 0. If this has not happened by the time k reaches M we 
increase M and repeat. 

For denominator 3 the value M = 25 is sufficient to show that 6(1/3) is finite for 3 < I < 2000. 
The following table shows what happens for the first few values. It gives the initial term, the 
number of steps to reach an integer, and the integer that is reached. 

t t .34 5 6789 10 U 

bldil - 33 3 3333 3 3 - '' 



reaches 



2 6 1 1 5 2 

1 8 1484710602474311520 2 7 8 3 1484710602474311520 220 
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(These are sequences A072340 and A085276 in [Sloane 1995-2003].) In the range I < 2000 large 
values of 9(1/3) are scarce. The first few record values are 0(1/3) = 0,2,6,22,23, reached at 
I = 3, 4, 5, 28, 1783 respectively. 

Starting values r = take longer to converge-we discussed the cases d < 8 in Section 1 . The 
initial values of 8((d + l)/d) can be found in sequence A073524 in [Sloane 1995-2003]. The first few 
record values are 0, 1, 2, 3, 18, 26, 56, 79, 200, 225, 388, 1444, reached at d = 1, 2, 3, 4, 5, 11, 19, 31, 37, 67, 
149, 199 respectively (sequences A073529, A073528). R. G. Wilson, v. [Wilson 2002] has checked 
that 6((d+ l)/d) is finite for d < 500. 

It is amusing to note that the record value 1444 has the following interpretation: starting with 
|SS and repeatedly approximately squaring, the first integer reached is 



a number with about 10 435 digits. 

The approximate multiplication map f r (x) of (43) is easier to compute since it grows more 
slowly. Let 9 T (n) denote the smallest value of k > 1 for which fr k \n) is an integer. We give just 
one example. This table shows what happens when f±/z(n) is iterated with starting value n: 

n: 1 2 3 4 5 6 7 8 9 10 11 12 ••■ 
6»4 /3 (n) : 13212918 3 1 7 2 1 ■■■ 
reaches : 4 4 4 8 84 8 84 20 12 84 20 16 

(sequences A085068 and A085071). Large values of #4/3(71) are again scarce. The first few record 
values are #4/3(71) = 1, 3, 9, 15, 17, 18, 24, 27, 28, 30, 40, reached at 0, 1, 5, 161, 1772, 3097, 3473, 23084, 
38752, 335165,491729 respectively (sequences A085328 and A085330). We thank J. Earls [Earls 
2003] for computing the last six terms in these two sequences. 

All the evidence supports the conjectures made here; it would be nice to know more. 
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